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GLOBAL WELL-POSEDNESS OF THE THREE-DIMENSIONAL VISCOUS 
PRIMITIVE EQUATIONS OF LARGE SCALE OCEAN AND ATMOSPHERE 

DYNAMICS 
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Abstract. In this paper we prove the global existence and uniqueness (regularity) of strong solutions 
to the three-dimensional viscous primitive equations, which model large scale ocean and atmosphere 
dynamics. 
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1. Introduction 

Large scale dynamics of oceans and atmosphere is governed by the primitive equations which are de¬ 
rived from the Navier Stokes equations, with rotation, coupled to thermodynamics and salinity diffusion- 
transport equations, which account for the buoyancy forces and stratification effects under the Boussinesq 
approximation. Moreover, and due to the shallowness of the oceans and the atmosphere, i.e., the depth of 
the fluid layer is very small in comparison to the radius of the earth, the vertical large scale motion in the 
oceans and the atmosphere is much smaller than the horizontal one, which in turn leads to modeling the 
vertical motion by the hydrostatic balance. As a result one obtains the system (1)—(4), which is known 
as the primitive equations for ocean and atmosphere dynamics (see, e.g., [20],[21],[22], [23], [24], [33] and 
references therein). We observe that in the case of ocean dynamics one has to add the diffusion-transport 
equation of the salinity to the system (l)-(4). We omitted it here in order to simplify our mathematical 
presentation. However, we emphasize that our results are equally valid when the salinity effects are taking 
into account. 

Let us remark that the horizontal motion can be further approximated by the geostrophic balance when 
the Rossby number (the ratio of the horizontal acceleration to the Coriolis force) is very small. By taking 
advantage of these assumptions and other geophysical considerations several intermediate models have 
been developed and used in numerical studies of weather prediction and long-time climate dynamics (see, 
e.g., [4], [7], [8], [22], [23], [25], [28], [29], [30], [31] and references therein). Some of these models have also 
been the subject of analytical mathematical study (see, e.g., [2], [3], [5], [6], [9], [11], [12], [13], [15], [16], 
[17], [26], [27], [33], [34] and references therein). 

In this paper we will focus on the 3 D primitive equations in a cylindrical domain 


n = M x (-h, 0), 
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where M is a smooth bounded domain in R 2 : 


dv dv -* 

—— \- (v • V)v + w— —b Vp + fk x v + Liv = 0 
dt dz 

(1) 

d z p + T = 0 

(2) 

V • v + d z w = 0 

(3) 

dT dT 

~KT+v ■ VT + w— + L 2 T = Q 

dt dz 

(4) 


where the horizontal velocity field v = ( V\,V 2 ), the three-dimensional velocity field (v\,V 2 ,w), the temper¬ 
ature T and the pressure p are the unknowns. / = fo(P + y) is the Coriolis parameter, Q is a given heat 
source. The viscosity and the heat diffusion operators L\ and L 2 are given by 


1 A 1 92 

Re 1 Re 2 dz 2 ’ 

1 1 d 2 

Rt\ Rt -2 dz 2 ’ 


(5) 

( 6 ) 


where Re\, Re 2 are positive constants representing the horizontal and vertical Reynolds numbers, respec¬ 
tively, and Rt\,Rt 2 are positive constants which stand for the horizontal and vertical heat diffusivity, 
respectively. We set V = ( d x ,d y ) to be the horizontal gradient operator and A = <9 2 + <9 2 to be the 

horizontal Laplacian. We observe that the above system is similar to the 3 D Boussinesq system with the 

equation of vertical motion is approximated by the hydrostatic balance. 

We partition the boundary of O into: 

T„ = {(x,y, z) G 0 : z = 0}, (7) 

T b = {(x,y, z) G O : 2 = -h}, (8) 

T s = {(x, y,z) GO : (x, y) G dM , —h < z < 0}. (9) 


We equip the system (l)-(4) with the following boundary conditions - with wind-driven on the top surface 
and non-slip and non-heat flux on the side walls and bottom (see, e.g., [20], [21], [22], [24], [25], [28],[29], 
[30]): 


r dv dT 

on T„ : — = hr, w = 0, — = -a(T - T )■ 
oz oz 

dv dT 

on r 6 : — = 0, u> = 0, — = 0; 

dz dz 

_ dv dT 

on T s : v ■ n = 0, — x n = 0, —— = 0, 
dn dn 


( 10 ) 

( 11 ) 

( 12 ) 


where r(x,y) is the wind stress on ocean surface, n is the normal vector to T s , and T*(x,y) is typical 
temperature distribution of the top surface of the ocean. For simplicity we assume here that r and T* are 
time independent. However, the results presented here are equally valid when these quantities are time 
dependent and satisfy certain bounds in space and time. 

Due to the boundary conditions (10)-(12), it is natural to assume that r and T* satisfy the compatibility 
boundary conditions: 


t • n = 0, 
dT* 


dr _ 

■kz; x n = 0 , 

dn 


dn 


= 0 on dM. 


on dM. 


(13) 

(14) 
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In addition, we supply the system with the initial condition: 

v(x,y,z,0) = v 0 (x,y, z). (15) 

T(x, y, z, 0) = Tq(x, y, z). (16) 

In [20], [21] and [33] the authors set up the mathematical framework to study the viscous primitive 
equations for the atmosphere and ocean circulation. Moreover, similar to the 3 D Navier-Stokes equations, 
they have shown the global existence of weak solutions, but the question of their uniqueness is still open. 
The short time existence and uniqueness of strong solutions to the viscous primitive equations model was 
established in [15] and [33]. In [16] the authors proved the global existence and uniqueness of strong 
solutions to the viscous primitive equations in thin domains for a large set of initial data whose size 
depends inversely on the thickness of the domain. In this paper we show the global existence, uniqueness 
and continuous dependence on initial data, i.e. global regularity and well-posedness, of the strong solutions 
to the 3 D viscous primitive equations model (1)-(16) in general cylindrical domain, O, and for any initial 
data. It is worth stressing that the ideas developed in this paper can equally apply to the primitive 
equations subject to other kids of boundary conditions. As in the case of 3 D Navier-Stokes equations the 
question of uniqueness of the weak solutions to this model is still open. 

2. Preliminaries 
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By integrating equation (2) we obtain 


P(x,y,z,t) = - / T(x,y,£,t)d4+p s (x,y,t). 

J-h 

Substitute (17) and the above relation into equation (1) we reach 

dv f f z \ dv 

— + (u • V)u - ( J V • v(x, y, £, t)dn — 

+Vp s {x,y,t)-X7 [ T(x,y,£,t)d£ + fk x v + L x v = 0. (25) 

J-h 

Remark 1. Notice that due to the compatibility boundary conditions (13) and (14) one can convert the 
boundary condition (10)-(12) to be homogeneous by replacing (v,T) by (v + ( z + h } ~ h / 3 r, T + T*) while 
(23) is still true. For simplicity and without loss generality we will assume that r = 0,T* = 0. However, 
we emphasize that our results are still valid for general r and T* provided they are smooth enough. In a 
forthcoming paper we will study the long-time dynamics and global attractors to the primitive equations 
with general r and T*. 


Therefore, under the assumption that r = 0, T* = 0, we have the following new formulation for system 

( 1 )—( 16 ): 


dv 

dt 


+ L x v + (v ■ V)v - (^J V • v(x, y, £, t)d^j 


dv 

dz 


+Vp s (x,y,t) -V / T(x,y,€,i)d£ +fkx v = 0, 

J-h 

dT 

— + L 2 T + vVT- 
dt 


f z \ dT 

J V • v(x,y,£,t)d£j — = Q, 


dv 

dz 


2=0 


= 0 , J 


z=—h 


dv _ 

= 0, v ■ n r =0, — x n 
dn 


= 0, 


(d z T + aT) | z=0 = 0; d z T\ z= _ h = 0; d n T | r> = 0, 
v(x,y,z, 0) = v 0 {x,y,z), 

T{x, y, z, 0) = T 0 (x,y,z). 


(26) 

(27) 

(28) 

(29) 

(30) 

(31) 


2.2. Properties of v and v. By taking the average of eciuations (26) in the z direction, over the interval 
{—h, 0), and using the boundary conditions (28), we obtain 


dv ( f z \ dv 

— + (n • V)n - ( J V • v(x, y,£, t)d£J — + Vp s (x,y,t ) - V 

+fk xv— —^—Av = 0 . 

Re i 


lj J T(x,y,£,t)d£dz 


As a result of (22), (23) and integration by parts we have 


(32) 


(v ■ V)i> 



V ■v(x,y,^,t)d4 


dv 

dz 


(■v ■ V)i> + [(i> • V)i> + (V • v) v\. 


(33) 
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By subtracting (32) from (26) and using (33) we get 


dv 

dt 

-V 


+ L\v + (v ■ V)u ~ (^J ^ ' v i x > 2/> £> t)d^j zp + {v ■ V)u + (v ■ V)i 


T(x,y,£,t)d£ - - 


T(x, y, £, t)dt;dz ) + fk x v = 0. 


— h J —h 


[(u • V)u + (V -v) v\ 


(34) 


Therefore, v satisfies the following equations and boundary conditions: 
dv 1 


dt Re i 


Av + (v ■ V)u + [(u • V)v + (V • v) v\ + fk x ' 


+V 


p a (x,J T ( x >y’*)^ dz 


V • v = 0, 


' —h J —h 
in M, 


= 0, 


dv 


v • n = 0, — x ft = 0, 

dn 


on dM, 


and v satisfies the following equations and boundary conditions: 


dv 

dt 


+ L\v + (v ■ V)v — 


h 


dv 

v • v(x, y , £, t)d£ ) -p + (v ■ V)v + (v ■ V)v 


-[(« • V)u + (V • v) v\ + fk x v - V ( j T(x,y,£,t)d4 - - J J 


T(x,y,£,t)d£dz ) = 0, 


dv 

dz 


2=0 


= 0, f 


z=—h 


~ dv _ 

= 0, v ■ n L = 0, — x n 
dn 


= 0 ,. 


(35) 

(36) 

(37) 


(38) 

(39) 


Remark 2. We recall that by virtue of the maximum principle one is able to show the global well-posedness 
of the 3 D viscous Burgers equations (see, for instance, [19] and references therein). Such an argument, 
however, is not valid for the 3 D Navier-Stokes equations because of the pressure term. Remarkably, the 
pressure term is absent from equation (38). This fact allows us to obtain a bound for the L 6 norm of v, 
which is a key estimate in our proof of the global regularity for the system (1)—(16). 


2.3. Functional spaces and Inequalities. Let us denote by L 2 (fl),L 2 (M) and H m (Q), H m (M) the 
usual L 2 —Lebesgue and Sobolev spaces, respectively ([1]). We denote by 


Let 


Vi 


V 2 



(In V> z )\ p dxdydz ) p , 
\ P dxdy)* , 


for every (j> G L p (fl) 
for every £ L P (M). 


(40) 


v £ C°°(n) : 


dv 

dz 


= 0 . A 

= 0 , 

n dv -* 
v ■ n L = 0, —— x n 

= 0. V ■ v = 0 ] 

O 

II 

z=—h 

lr * dn 

r. J 


T £ C°°{n) : 


dT 

dz 


z=—h 


(dT 

~° ; \-d~z +aT 


= 0 ; 

2=0 



We denote by Vi and V 2 be the closure spaces of Vi in H 1 ( fi), and V 2 in under H 1 — topology, 

respectively. 
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Definition 1. Let Vo G V\ and To G Vi, and let T be a fixed positive time, (v, T) is called a strong solution 
o/(26)-(31) on the time interval [0,7] if it satisfies (26) and (27) in weak sense, and also 

vGC([0,T],Vi)nL\[0,T],H 2 (n)), 

T G C([0, T\, V 2 ) n L 2 ([ 0 , T], TT 2 (fi)), 

^GL\[0,T],L 2 m, 

^ GL\[0,T],L 2 m). 

For convenience, we recall the following Sobolev and Ladyzhenskaya’s inequalities in R 2 (see, e.g., [1], 
[10], [14], [18]) 

j^||L4(M)<Q,||<(>||^ 2 ||0||^ M) , (41) 

ll^llLS(M) < C'oII^II^mjII^IIhmm)’ ( 42 ) 

for every (f £ H 1 (M), and the following Sobolev and Ladyzhenskaya’s inequalities in R 3 (see, e.g., [1], [10], 

[14], [18]) 

IIV’IMfi) < C' 0 ||w||^ / 2 2 (n) ||u||^ / 1 2 (n) , (43) 

l u llL 6 (n) < ColMI/iHn), (44) 

for every u £ TT 1 (fl). Here Co is a positive constant which might depend on the shape of M and H but 
not on their size. Moreover, by (41) we get 

IHIIV) = IIH 3 || 4 L*(M) < Co|||0| 3 ||| 2( M)lll^| 3 || 2 ffl (M) 

< C-o||0||i« (M) \fi\ A |V</>| 2 dxdy^j + \m 12 e (M) , (45) 

for every (f> £ TT 1 (Af). Also, we recall the integral version of Minkowsky inequality for the L p spaces, p > 1. 
Let Hi C R mi and H 2 C R m2 be two measurable sets, where m\ and m 2 are two positive integers. Suppose 
that is measurable over Hi x H 2 . Then, 



\f(t,v)\dv] df 


1 1 /p 


< 


(jf \m,ri)\ p d£) /P dy. 


(46) 


3. A Priori estimates 

In the previous subsections we have reformulated the system (1)—(16) and obtained the system (26)-(31). 
The two systems are equivalent when (v,T) is a strong solution. The existence of such a strong solution 
for a short interval of time, whose length depends on the initial data and the other physical parameters of 
the system (1)—(16), was established in [15] and [33]. Let (vo,T 0 ) be a given initial data. In this section we 
will consider the strong solution that corresponds to this initial data in its maximal interval of existence 
[0,7),). Specifically, we will establish a priori upper estimates for various norms of this solution in the 
interval [0,7)). In particular, we will show that if 7) < 00 then the H 1 norm of the strong solution is 
bounded over the interval [0,7)). This key observation plays a major role in the proof of global regularity 
of strong solutions to the system (1)—(16). 
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3.1. L 2 estimates. We take the inner product of equation (27) with T, in L 2 (U), and obtain 

l^r + i^ l|VT|1 ^ + w 2 ml + amz = 0)l1 ^ 


= / QT dxdydz — I ( v • VT — 

Jo, Jo, 


/ z \ \ 

V • v(x, y, £, t)d£] —) T dxdydz. 


After integrating by parts we get 

I ^\J ^ ' v ( x, y’ -g-^J T dxdydz = 0. 

As a result of the above we conclude 

+ w} vnl + w} T ^ 1 + amz = 0)1,2 

= f QT dxdydz < ||Q|| 2 ||T|| 2 . 

J n 

Notice that 

\\T\\l<2h 2 \\T z \\l + 2h\\T(z = 0)\\l 
Using (48) and the Cauchy-Schwarz inequality we obtain 


+ w^ nl + w 2 nl+amz = 0)lli 


<2 (h 2 Rt2 + -)\\Q\\l 

a 

By the inequality (48) and thanks to Gronwall inequality the above gives 

\\T\\l < e- ||To||i + (2 h 2 Rt 2 + 2h/a) 2 \\Q\\ 2 , 

Moreover, we have 

f [t^I|VT( s )||2 + J_ \\T z (a)\\l + all T{z = 0)(«)|||1 ds 
J o L-n3l -txi 2 

< 2 {h 2 Rt 2 + -)\\Q\\l t + e~ ||Xb||l + (2 h 2 Rt 2 + 2V«) 2 1|Q || ^ 

a 

By taking the inner product of equation (26) with v, in L 2 (0), we reach 


1 d||u|| 


1 + 


2 dt Re i 

= -[ {v-V)v-(f V • v(x, y, £, t)d£ 
Jn . \J—h 


dv 

dz 


■ v dxdydz 


+ J yfk x v + Vp s — V yj T(x,y, £, t)d£ ) ) • v dxdydz. 


By integration by parts we get 

[ (v’V)v-(f V -v(x,y,£,t)d4 
■In . \J—h 


dv 

dz 


■ v dxdydz = 0. 


/ Vp s • v dxdydz = h Vp s • v dxdy = —h / p s (V • v) dxdy = 0. 

Jn J m J n 


By (36) we have 
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Since 


(fk x v) ■ v = 0, 


then from (53)-(55) we have 


1 d\\v 


2 dt 


Re\ 


IIVwlll 


Re 2 ' 


u z 112 


in J-h 

<fcimi 2 iiv«ii 2 . 

By Cauchy-Schwarz and (51) we obtain 


f f T(x,y,£,t) d£(V • v) dxdydz 

Jn J-h 


ld|H| 


1 


1 




2 dt Re i" " z Re 2 ' 

< h 2 Re 1 ||T|| 2 < h 2 Re x (||T 0 || 2 + (2 h 2 Rt 2 + 2h/a) 2 \\Q\\ 2 2 ) . 

Recall that (cf., e.g., [14] Vol. I p. 55) 

\\v\\l<C M \\S7v\\ 2 . 

By the above and thanks to Gronwall inequality we get 

IMl! < e “ c ' M " ei (*IIW||l + 11 ^ 0 1| 2 ) 

+C M h 2 Re\ [ ||Tb||l + (2 h 2 RU_ + 2h/a) 2 \\Q\\ 2 } . 

Moreover, 


-^-||Vw(s)||| + | 2 


ds 


< h 2 Re\ (||To||l + (2 h 2 Rt 2 + 2*./«) 2 1|Q|| 1) t + e“ 
+C M /* 2 i?e 2 [ 11To111 + (2 h 2 Rt 2 + 2V«) 2 ||Q||1 ] • 
Therefore, by (51), (52), (56) and (57) we have 


mil 


+ \\T(t)\\l 


Re 1 


l|Vn(s)||( 


R&2 


ds 


^l|VT( S )||2 + -^||T z ( s )||2 + a||T(z = 0)( S )|||| ds < K^t), (58) 


( 55 ) 


(56) 


(57) 


where 


Ki(t) = 2(h 2 Rt 2 + h/a)\\Q\\l;t+ (h\\v 0 \\l + ||n 0 ||i) 

+ (l + Cm h 2 Re 2 + h 2 Re 1 1) [ ||To || 2 4" (2 h 2 Rt 2 + 2hja) 2 1| Q || \ ] . (59) 

3.2. L 6 estimates. Taking the inner product of the equation (38) with |n| 4 n in L 2 (fl), we get 


ld||i 


1 


6M + -ik 1 ,|wh ” 


J (jVn| 2 |n| 4 + |VIXJ1 2 1“ |u| 2 ) dxdydz + — j (jt>z| 2 |v z | 4 + |d z |n| 2 | 2 |n| 2 ^ dxdydz 


(v ■ V)i> — 
+fk x v — V 


dv 


V • v(x, y, £, t)dS, ) ——{- (v ■ V)i> + (v ■ V)v — [(i> • V)v + (V • v) v\ 


— h 


T(x, y , £, t)d£ - - j I T(x, y, f, t)d^dz ) \ • \vyv dxdydz. 


— h J —h 
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By integration by parts we get 


• |r)| 4 f dxdydz = 0. 


(v ■ V)i> — 


dv 

1 V • v(x,y,£,t)d£j -rp 


Since 


(/fc x v) 


\v\ v = 0, 


then by (36) and the boundary condition (28) we also have 


/ (v • V)v • |u| 4 i> dxdydz = 0. 

Jn 


Thus, by (60)-(62) we have 


+ _L ^ ^|Vu| 2 |u| 4 + |VM 2 |" |tf) dxdydz + (\v z \ 2 \v z \ 4 + |<9>| 2 |~ [i 

= — | (v ■ V)rJ — (v • V)u + (V • u) v 


-V { / T(x,y,£,t)d(, — — 


J j T(x,y,£, t)d^dzj | • \v\ 4 v dxdydz. 


Notice that by integration by parts and boundary condition (28) we have 


/ n L 


(v ■ V)v — [(u • V)v + (V • v) v\ 

•0 pz 


-V ( f T(x,y,£,t)d£ ~ t [ [ T(x,y,£,t)d£dz 

\J —h n J-hJ-h 


\v\ 4 v dxdydz 


In L 


(V • v) v ■ |u| 4 w + (v • V)(|t;| 4 u) • v — v k vi d Xk (\v\ 4 v ■ 7 ) 

0 nz 

-h 10 J —h J — h 


T(x, y, £, t)d£ - ^ J ^T(x,y,£, t)d£dz \ V • (|n| 4 u) 


dxdydz. 


(60) 


(61) 


(62) 


| 2 ^j dxdydz 
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Therefore, by Cauchy-Schwarz inequality and Holder inequality we obtain 

MHle . 1 f / |Y7 ~ |2 |~|4 , |x-7|~|2| 2 I~|2\ , 1 [ ( |~|2i~ 


J ^|Vu| 2 |u| 4 + |V|u| 2 |" |u| 2 ^ dxdydz + J (\v z \ 2 \v z \ 4 + |9 z |t)| 2 | 2 |v| 2 ) dxdydz 


< C / |v| / |Vi>| M 5 dz dxdy 

JM L J-h 


|Vt>| |u| 4 dz ) dxdy 


JM L \J-h / \J-h / J 

+C f W\ f |Vw| |u| 4 dz dxdy 

JM [ J-h 

< C J 1,| (/ |Vu| 2 \v\ 4 dz^j (/ \v\ 6 dz'j dxdy 


v\ 2 dz) [ / |Vv| 2 |u| 4 dz ) ( / |u| 4 <ij) dxdy 


JM l\J-h J \J-h J \J-h J 

+C [ [|T| f f |Vu| 2 |u| 4 d,^ ( [ |w| 4 dz] dxdy 


< C\\v\\ L i(M) ^ |Vu| 2 |u| 4 dxdydz ^ (^j 


u| 6 dz dxdy 


4 \ i/4 


u| 2 dz dxdy 


| Vd| 2 M 4 dxdydz 


2 \ l/ 4 


v | 4 dz dxdy 


+C\\ \T\ \\ LHM) |Vu| |u| 4 dxdydz J 
By using Minkowsky inequality (46), we get 


o \ 2 \ 1/4 

|u| 4 dz ) dxdy J 


2 \ 

u| 6 dz ) dxdy I < C 


|u| 12 dxdy dz 


By (45), 


f |u| 12 dxdy < Co ( I |u| 6 dxdy] f f |u| 4 |Vu| 2 dxdy \ + ( [ M 6 dxdy] . 

JM \Jm ) \Jm J \Jm J 


Thus, by Cauchy-Schwarz inequality we obtain 

/ „ , „n \ 2 \ 1/2 


dz^j dxdy j < C\\v\\ 3 L6(n) (^J^v\ 4 \Vv\ 2 dxdydzj + |M|£ 6(n) . (63) 


Similarly, by (46) and (42), we also get 

(nr |u| 4 dz^ dxdy] < C 


\ I/ 2 

|i>| 8 dxdy I dz 


< C / N||6 (m) (||Vu|| i2(M ) + \\v\\l H m)) dz < C||u||l (||W|| 2 + ||u|| 2 ), 
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0 \ 4 

IJTI2 \ 


vf dz I dxdy < C 


II L 6 (M) + IMI L 2 (M) 


I?;| 8 dxdy dz 


) dz < C||S||= /2 (||V»|| 


i' 2 + IlSIlf 


Therefore, by (63)-(65) and (41), we reach 

(l V ^ 2 |^| 4 + l V l^| 2 | 2 l^l 2 ) dxdydz + P^J Q (l^| 2 |^| 4 + |^N 2 f |v| 2 ) dxdydz 
. r \ 3 / 4 


< C\\v\\\ /2 ||Vn||2 /2 ||c||6 /2 ( / |Vu| 2 M 4 dxdydz ) + C\\v\\-^ ||VuH^^Hvlle 




+C1MII (l|Vv|| 2 + IM| 2 ) / |Vu| 2 |v| 4 dxdydz 


+C\\T\\l /2 ||VT||2 /2 ||u||g /2 (||Vz;||2 /2 + |H|2 /2 ) (J^ IVu| 2 |£| 4 dxdydz 
Thanks to the Young’s and the Cauchy-Schwarz inequalities we have 

^IMIl , 1 [ /|VT~|2|~|4 , I V7l~l2l 2 |~|2\ j i. , 1 [ ( |~|2|~ 


dt Re i J n V 

• /^Il7rrl|2 llT7T7l|2||~||6 


J (|Vu| 2 |?;| 4 + |v|u| 2 |" H 2 ) dxdydz+ rt— j (\v z \ 2 \v z \ 4 + \d z \v\ 2 \ 2 \v\ 2 ^ dxdydz 


n N 

l|2||~l|6 


< C\\v\\ 2 || Vu|| 2 ||u|| 8 + C||u||i||W|| 2 + C||T|| 2 ||VT|| 2 + C||u|| 2 ||u|||. 
By (58) and Gronwall inequality, we get 


/ ((Re7 J l Vu | 2 |^| 4 dxd y dz + J l^l 2 !^! 4 dxdydz ^ < K 6 (t), 


where 


K 6 (t) = e K ? w 


Taking the inner product of the equation (27) with |T| 4 T in L 2 (fl), and by (27), we get 

+ 4rJ n ,VT,2,T i 4 dxdydz + w 2 j n |Tz|2|T|4 dxdydz +a||T(z = 0)1,8 

= J Q\T\ 4 T dxdydz — J • X7T — "V ■ v(x,y,^,t)dpj \T\ 4 T dxdydz. 
By integration by parts and (36) we get 

-J (v-VT-^J V ' v ( x > V> £> tfdPj \T\ 4 T dxdydz = 0. 

As a result of the above we conclude 


Mmn 

6 dt 


+ 4r I NT\ 2 \T\ 4 dxdydz +[ \T Z \ 2 \T\ 4 dxdydz+ a\\T(z = 0)\\ 6 6 
■Kti Jn tit2 Jn 


= / Q\T\ 4 T dxdydz < \\Q\\ 6 \\T\\l 


By Gronwall, again, we get 


||7V)I| 6 < ||Q||ffi(f2) t+ ||T 0 || ff i(n). 
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3.3. H 1 estimates. 

3.3.1. 11VU11 2 estimates. First, let us observe that since v is a strong solution on the interval [0,7),) then 
Av € L 2 ([0, %), L 2 (M)). Consequently, and by virtue of (36), Av ■ n € L 2 ([0,7)), if ” 1/,2 (<9M)) (see, 
e.g., [10], [32]). Moreover, and thanks to (36) and (37), we have Av ■ n = 0 on dM (see, e.g., [35]). This 
observation implies also that the Stokes operator in the domain M, subject to the boundary conditions (37), 
is equal to the —A operator. 

As a result of the above and (36) we apply a generalized version of the Stokes theorem (see, e.g., [10], 
[32]) to conclude: 


/ Vp s (x, y, t) ■ Av(x, y, t)dxdy = 0. 

J M 


By taking the inner product of equation (35) with —Av in L 2 (M), and applying (36) and the above, 
we reach 


1 d||W|| 

2 dt 


+ -^—]|An||i = f {(n • V)v + [(n • V)n + (V • v) n]} • An dxdy + [ fkxv-Avdxdy. 
fiei ./ i\/r f ' .1 


Re\ j M >- J J m 

Following similar steps as in the proof of 2D Navier Stokes equations (cf. e.g., [10], [32]) one obtains 

<C||U||J /2 ||Vn|| 2 ||An||^ /2 . 


/ ( v ■ V)v ■ Av dxdy 

JM 


IM 

Applying the Cauchy-Schwarz and Holder inequalities, we get 

,o 


f (v ■ V)v + (V • v) v ■ Av dxdy < C f f |zi| |Vu| dz |Az;| dxdy 

JM ' J M J — h 

dxdy 


< C 


< C 


IM 


IM J-h 

r 0 \ 1 / 2 / f -0 \ !/ 2 

|n| 2 |VF| dz) ( |Vn| dz ) |Au| 


-h 
•0 

M \J-h 


- 1/4 

r / p0 \ ^ 

1/4 

r 


/ ( / |Vv| dz ) dxdy 


/ Av 2 dxdy 


JM \J — h ) 


J M 


n 1/2 


< C , ||VF|| 2 /2 (^j |n| 4 |VF| 2 dxdydz^j ' ||AtJ|| 2 . 

Thus, by Young’s and Cauchy-Schwarz inequalities, we have 

d ^ d ^ 2 + + C'llVulla + C J n l^| 4 l W l 2 dxdydz + C\\v\\ 2 2 . 

By (58), (66) and thanks to Gronwall inequality and we obtain 


\\Vv\\ 2 2 + \Av\ 2 2 ds < 


(70) 


where 


A 2 (t) — e Kl ^ ||uo||^i(fj) + Ki(t) + Ke(t) 


( 71 ) 
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3.3.2. |H| 2 estimates. Denote by u = v z . It is clear that u satisfies 

du / f z \ 9a 

— + Lm + (v ■ V)u - ( J V • v(x, y , £, t)d£J — 

+(u • V)i> — (V • v)u + fk x u — VT = 0. (72) 

Taking the inner product of the equation (72) with u in L 2 and using the boundary condition (28), we get 

= — J ^(i> • V)u — (/ V • v(x, y, £, t)dt^J — ^ • u dxdydz 
— J ^( u ■ V)u — (V • v)u + fk x u — VT^j • u dxdydz. 

By integration by parts we get 

— J ^(v-V)u-^J V • v(x, y, £, t)dt^J • u dxdydz = 0. (73) 


(fk x u) ■ u = 0, 


then by (73) and (74) we have 


= — ((u ■ V)v — (V • v)u — VT) • u dxdydz 

Jn 

< C J (|v|) |m| |Vu| dxdydz + ||T|| 2 ||Vu|| 2 

< C'lK'lle ||«|| 3 ||v«|| 2 + ||T|| 2 ||Vu|| 2 
<C|H| 6 ||«|| 2 /2 ||V«||* /2 + ||T|| 2 ||Vu|| 2 . 

By Young’s inequality and Cauchy-Schwarz inequality, we have 


<C1Mle ||w|| 2 + C||T||2 
<C'(||Vu||^ + rt) \\u\\l + C\\T\\l. 


By (58), (66), (70), and Gronwall inequality, we get 


where 


\\vz\\l + ~jT f l|V^(a)||l + -±-f ||t7*,(a)||l ds < K z (t), 
Rei Jo Re 2 Jo 

K z (t) = e(*2M+*. a/8 W>* [||« b ||| rl(n) + K,(t) . 
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3.3.3. 11 Vz'11 2 estimates. By taking the inner product of equation (26) with —An in L 2 (fl), we reach 


IdJMf 1 a J_ 112 

2 dt + V " 2 W *" 2 


t)d£ 


[ (v • V)u - ( [ V • v(x,y,£, 
Jn L \J-h 


+fk x n + — V ( / T(x,y,£,t)d£ 


dv 

dz 


i-h 


Av dxdydz 


p ‘ pO pO 

<C \v\\Vv\ + |Vu| dz\v z \ + / |VT| dz 
Jfl . J —h J—h 


| At) | dxdydz 


<C\\v\\ Le(n) \\Vv\\ L 3 {n) \\Av\\ 2 + C / \Vv\dz \v z \\Av\ dz ) dxdy + C\\VT\\ 2 \\Av\\ 2 . 

Jm \J —h J-h 


Notice that by applying the Proposition 2.2 in [5] with u = v, f = Av and g = v z , we get 


[ (f \Vv\dz ( \v z \\Av\ dz) dxdy <C\\Vv\\l /2 \\v z \\l /2 \\Vv z \\l /2 \\Av\\l /2 . 
Jm \J—h J-h J 


As a result and by (43) and (44), we obtain 


l^^ + itJ Avlll + it 2 llVVz ^ 


< 


c (|MUe(Q) + ||Vu|| 2/2 K|| 2/2 ) ||Vu|| 2/2 ||Au||^ /2 + /i||VT|| 2 ||Au|| 2 . 


Thus, by Young’s inequality and Cauchy-Schwarz inequality, we have 




< 


c(|Hlle (n) + ||Vn|| 2 |K|| 2 ) ||Vn|| 2 + C||VT|| 


By (58), (66), (70), (75) and thanks to Gronwall ineciuality, we obtain 


||Vu|| 2 +/ (^||A«( 8 )||£ 


Re 2 


Vz;^(s)||| ) ds < Ky{t), 


where 


K v (t) = e^6 /3(t) t+K RV K R t '> 


m( o.) 


■ K\ (t) 
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3.3.4. ||T||h-i estimates. Taking the inner product of the equation (27) with —AT — T zz in L 2 ( fl), we get 
ld(||VT|||+||T z ||| + a||VT( 2 = 0)||I) 


dt 


+^I|VT 2 ||2 + (^ + ^) (l|VT,||i+a||Vr(« = 0)||i) + -^-||T s 


2 

zz ||2 


/a L 


v • VT- 


'-h 


V • « df T 2 - Q 


[AT + T zz ] dxdydz 
r° 

l-h 


<C [ (M|VT| + |Q|) |AT + T 22 | dxdydz + [ \f \Vv\ dz [ \T Z \ |AT + T 22 | dz 
Jn JM IJ-h J—h 


dxdy 


1/2 


<C|H| 6 IIvt || 3 (||at||| + ||vr 2 ||| + ||t^||) 

+C||V w ||^ /2 ||Ai;||^ 2 ||T 2 ||^ /2 (||Ax’ll! + IIVX.II! + ||T 22 ||!) 3/2 + ||Q ||2 (||AX|| 
< C [|M| 6 l|VX||V 2 + UVull^llAell^llX.ll^l (||AT|| 2 + ||VX.||! + ||T 22 |||) 
+ ||Q|| 2 (||AT||2 + ||VT 2 ||2 + ||T 22 ||2) 1/2 . 


|vt 2 


\\Tzzwi) 


211/2 


\ 3/2 


By Young’s inequality and Cauchy-Schwarz inequality we have 

d(l|Vr||l + ||r z ||l + q||VX(^ = 0)||1) 

dt 

+ ^I|VT 2 ||2 + (||VX 2 ||! + a||VT(s = 0)|||) + ^||T„||1 

< C (Ml + 1|Vw||!iiAu||!) (11vx11! + ||t 2 ||!) + q|Q|||. 


By (66), (77), and Gronwall inequality, we get 
12 ' 11X,11! + a||VX(z = 0)||! 


J) VX|| 2 
r t r 




VT 2 


\l + 


1 

Rti 


1 

Rt.2 


I vx 


■a||VX(z = 0)||!)— IIT 


2 

zz \\2 


ds < Kt, 


where 


t+Ky(t) 



2 

m(n) 



(79) 

(80) 


4. Existence and Uniqueness of the Strong Solutions 

In this section we will use the a priori estimates (58)-(79) to show the global existence and uniqueness, 
i.e. global regularity, of strong solutions to the system (26)-(31). 

Theorem 2. Let Q £ H 1 ^), vo € V\, To G V 2 and T > 0, be given. Then there exists a unique strong 
solution (v,T) of the system (26)-(31) on the interval [0,T] which depends continuously on the initial data. 


Proof. As we have indicated earlier the short time existence of the strong solution was established in [15] 
and [33]. Let (v,T) be the strong solution corresponding to the initial data (vq,Tq) with maximal interval 
of existence [0, T„). If we assume that % < 00 then it is clear that 

lim sup (IMI/iqn) + ||r||tfi(fi)) = 00 . 

Otherwise, the solution can be extended beyond the time However, the above contradicts the a priori 
estimates (75), (77) and (79). Therefore % = 00 , and the solution (v,T) exists globally in time. 
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Next, we show the continuous dependence on the initial data and the the uniqueness of the strong 
solutions. Let (m,Ti) and ( 112 , 12 ) be two strong solutions of the system (26)-(31) with corresponding 
pressures (p s )i and (p s ) 2 , and initial data ((uo)lr(^o)i) and (( 110 ) 2 ,( 20 ) 2 ), respectively. Denote by u = 
Hi — ii 2 , q s = (p s ) 1 ^ (ps )2 and 8 = T\ — T 2 . It is clear that 


du 

Ut 


+ L\u + (111 • V)u + ( u • V)ii 2 - V • Hi (x,y, £, t)d%j ^ - (^j V • u(x, y, £, <)d^ 


<9u 2 

<9z 


+fk xu + Vq s -V 

BO 

— + L 2 8 + V1-S70 + U- VT 2 
at 


0 (%,y,t,t)d 4 = 0, 


J h V-Hi (x,y,£,t)d£j - (/ -u(x,y,^,t)d^j 


m 

dz 


= 0. 


u(x,y,z,t) = (n 0 )i - (h 0 ) 2 , 
6»(x, y, z, 0) = (T 0 )i - (T 0 ) 2 . 


(81) 


(82) 

(83) 

(84) 


By taking the inner product of equation (81) with with u in L 2 (fl), and equation (82) with 8 , in L 2 (fl) we 
get 


MHII 

2 dt 


Re 1 


IIV11II2 


i?e 2 


As || 2 


J (m • V)u + (u • V)h 2 - ^ V • v 1 (x,y,t,t)d£j ^ ^ V • u(x,y,£,t)d(^J ^ • u dxdydz 


In L 


fk x u + X7q s — V 


0{x,y,£,t)d€ 


'-h 


u dxdydz , 


and 


+ i ^ ||V0| > 2 + ^11**112+ a ii 0 ^ = °)iii 


^ v 1 -V 8 + u-VT 2 - (^J V • vi(x,y,£,t)d£j ^ V • u{x,y,£,t)d£ 


m 

dz 


By integration by parts, and the boundary conditions (28) and (29), we get 

du 
dz 


■u dxdydz = 0, 


[ ((v%-V)u-([ V • vi(x,y,£,t)d(, 

Jn \ \J—h v 

J V • v 1 (x,y,£,t)d£j ■ 8 dxdydz = 0. 


Since 


(/* 


A: x -u ) • = 0, 


0 dxdydz. 

(85) 

( 86 ) 

( 87 ) 


Then by (85), (86) and (87) we have 

1 


MMl 2 

2 dt 


+ 


Re 


■||Vu||( 


i?e 2 ' 


■**112 


= — f (u ■ V)n 2 • u dxdydz + I I V ■ u(x, y, £, t)d£,- 7 ^ • u dxdydz. 
Jn ' 2 n ./-ft, 02 
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l^r + + ik m * + aP{z = 0)1,2 


Notice that 


= — f (u ■ V)T 2 0 dxdydz + f [ V • u{x,y,^,t)d£^-^-6 dxdydz. 
Jn ' Jo. J-h oz 


(u ■ V)u 2 • u dxdydz < 11 Vz^ 2 11 2 11«11 3 117^ 11 6 < C || Vv 2 1| 2 11 tz || 2 /2 11 Vzz || , 

(u ■ W)T 2 d dxdydz < 11 Vi> 2 11 2 11$11 3 11ii11 6 < C||VT 2 || 2 ||0||^ /2 ||V0||^ /2 ||V U || 2 . 


Moreover, 


I / V • u(x,y,£,t)d!;—— ■ u dxdydz < / 
n J-h dz ./a 

I (r w-j f mA i!2 (f 


^ ( f \Vu\dz f \d z v 2 \ |u| dz\ dxdy 
M \J-h J-h ) 

0 xl/2\ 

|m| 02 ) dxdy 


\\/u\dz dxdy 


\d z v 2 \ 2 dz dxdy 


u\ 2 dz dxdy 


By Cauchy-Schwarz inequality, we get 


\S7u\dz dxdy < C||Vu|| 2 . 


By using Minkowsky inequality (46) and (41), we obtain 


u\ 2 dz dxdy < C 


|w | 4 dxdy dz 


<C f \u\\Vu\dz<C\\u\\ 2 \\Vu\\ 2 , 

J-h 


|<9 2 u 2 | 2 dz dxdy < C 


\d z v 2 \ A dxdy dz 


< C / \d z v 2 \\\/d z v 2 \ dz < c\\d z v 2 \\ 2 \\wd z v 2 \\ 2 . 

J-h 


Similarly, we have 


I / V • u(x, y, £, t)d£—-^-9 dxdydz < C\\ V M || 2 ||a 2 T 2 || 4/2 || Va 2 T 2 || 4 / 2 ||0|| 4/2 || Wd\\\ /2 . 

nJ-h dz 
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Therefore, by estimates (88)-(93), we reach 


i rf(Nll + l|fl|||) 
2 


dt 


+ 11 Vti 11 2 + -5-Klli + ^-||V0||1 + — ||0 S ||1 + all e(z = 0)111 


Re i 


Re 2 


Rt i 


Rt 2 ' 


< 


c ( || V ^2 1| 2 + ||^ 2 || 2 /2 ||Vd z u 2 ||' /2 ) |M| 1/2 ||Vu|] 


3/2 


+C'||VT 2 || 2 ||0||^ /2 ||V0||^ /2 ||V U || 2 + C||V U || 2 ||5 z T 2 ||^ /2 ||Va 2 T 2 ||^ 2 ||0||^ /2 ||V0||^ /2 . 
By Young’s inequality, we get 
d\\u n2 


dt 


< c(\\vv 2 \\i + 1|vTbiil + l|9 2 u 2 || 2 ||va z u 2 || 2 + ||a z T 2 || 2 ||va z T 2 ||l) (h 


0\\l 


Thanks to Gronwall inequality, we obtain 


ll«(*)lli + IIWII 2 <(IK* = o)||l + ll^ = o)|| 2 )x 

exp | cj* (|| Vu 2 (s)||| + ||VT 2 (s )||2 + ||d z u 2 ( S )|| 2 || Vd z u 2 ( S )|| 2 + ||3 Z T 2 ( S )|| 2 || va z r 2 ( s )|| 2 ) dsj . 
Since (v 2 ,T 2 ) is a strong solution, we have 

ll«(*)lli + II^WIll < (IK* = 0)111 + ||0(t = 0)111) exp{C (Kit + Kit + K Z K V + A) 2 )}. 


The above inequality proves the continuous dependence of the solutions on the initial data, and in partic¬ 
ular, when u(t = 0) = 0(t = 0) = 0, we have u(t) = 0(t) = 0, for all t > 0. Therefore, the strong solution 
is unique. 
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